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1. INTRODUCTION 


A need presently exists for better low-frequency noise 
reduction techniques. The f undamentai methods for noise 
attenuation such as transmission Joss or absorption are 
good for high frequency noise but provide limited success 
at reducing low-frequency noise. Analytical approaches to 
noise control also prove more effective for evaluating high 
frequency noise than low frequency noise. In order to 
better evaluate and attenuate low frequency noise, new 
analysis techniques are needed. 

A low frequency noise problem that often confronts 
noise control engineers in the transportation industries is 
cavity 'boom.' Cavity 'boom' is characteristic of an acous- 
tic resonance condition. The resonance occurs primarily as 
the result of acousto -structural coupling between the cav- 
ity and a compliant portion of the enclosing structure. 

The structure by its motion excites one or more of the low 
order acoustic modes. The source of the excitation is 
inherently difficult to identify and evaluate. 

Although the frequencies of the resonant modes in an 
acoustic cavity are determined by spectral analysis tech- 
niques, the character of the acoustic mode shapes is not 
currently obtained experimentally. For example, the two- 
microphone intensity method is ineffective for measurement 

i 

in reactive fields. A knowledge of the mode shapes, in 
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particular in terms of the particle motion, is advantageous 
in selecting an optimum noise attenuation solution. 

The primary objective for the research prepented here 
is the development of a procedure whereby the acoustic 
modal parameters, natural frequencies and mode shapes, in 
the cavities of transportation vehicles are determined 
experimentally. Furthermore, it should be possible to 
describe the acoustic mode shapes in terms of the particle 
motion. 

The well-known structural modal anlaysis techniques 
have long been adapted to mini-computer based dynamic 
analysis systems. A secondary objective for this study is 
to tailor the acoustic modal analysis procedure to existing 
mini-computer based spectral analysis systems. 
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2. ACOUSTIC THEORY 


2 . 1 Dynamic System Identification 

2 . 1.1 ZX£SM&n£X-E93S2n3JiL 

Harmonic excitation ie often encountered in engineer- 
ing systems. Although pure harmonic excitation is less 
likely to occur than periodic or other types of excitation, 
understanding the behavior of a system undergoing harmonic 
excitation is essential in order to comprehend how the sys- 
tem will respond to more general types of excitation. Har- 
monic excitation may be in the form of a force or displace- 
ment at some point in the system. 


The frequency response method is a harmonic analysis. 

A sinusoidal excitation is applied to a system and its 
steady-state response is examined over a frequency range of 
interest. For a linear system, both the excitation and the 
system response are sinusoidal of the same frequency. For 
a simple example, consider a damped, single degree-of- 
freedom system subject to harmonic excitation, Figure 1. 


\\ \\ \\ \\ 



Figure 1. Single DOF System with Forcing Function 
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The second-order linear differential equation of motion is 
written as [1]. 

mx (t) + cx(t) + kx(t) - f(t> 2.1 

where the forcing function is given by 

f(t) - F o e^ tJfc 2.2 

The excitation frequency, w, is sometimes referred to 
as the driving frequency. By substituting equation 2.2 
into equation 2.1 and dividing by m, the equation of motion 
becomes 

x (t) + 2Cw n x(t) + «*x(t). = ± F 0 e jWt 2.3 

where 


u 


n - ^ j~ - undamped natural frequency 


2.4 


and 


c 

C * 2maT* = dam Pi n 9 factor 2.5 

If a particular solution to equation 2.3 is assumed to be 
of the form 

x ( t ) * Xe^ 2.6 

an expression for the response is found to be 


V* 


l-(w/w n ) + j 2Cw/w n 


Xe 




2.7 


where 


ORIGINAL PA€i L:'j 
OF POOR QUALITY. 
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X - 


V* 


Ul-(G./w n )V + (2{«/« n ) 2 ) 1/2 


2.8 


and 


<P m tan 


“1 ZC<J/(i) n 
l-(w/w n ) ; 


2.9 


represent displacement amplitude and phase angle. 


The excitation is given by the real part of f(t) [2]. 
Then, the response is given by the real part of x(t). 
Therefore, the response will be regarded as the real part 
of the complex quantity satisfying equation 2.3. Thus, the 
steady-state response is shown to be 


x(t) 


Re 


F o e 


jwt 


A 


l-(o)/a» n ) 4 + j2fw/o> nJ 


2.10 


By evaluation of equation 2.10, it is seen that the 
response, x(t), is proportional to the force, f(t). The 
proportionality factor expressed in the frequency domain is 
given by 


H(w) 


1/k 


l-(w/w n ) 2 +j2Cw/w n 


2.11 


Equation 2.11 is known as the complex frequency response. 
The frequency response represents the ratio of the response 
to the excitation under steady-state conditions. Frequency 
response method is a mathematical mechanism for character- 
izing a dynamic system. It is noteworthy from equation 
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2.11 that if the driving frequency is near the natural fre- 
quency, the magnitude of the frequency response is limited 
only by the damping factor. 

Detailed analysis of the frequency response for a 
damped, single degree-of -freedom system appears in the text 
and will not be addressed here [3]. 

2.1.2 Transfer Function 

A fundamental concept employed by engineers for the 
purpose of characterizing a dynamic, linear system is that 
of system identification. As represented by the 'black 
box' analogy, the engineer by subjecting a system to a 
known input and recording the response can characterize the 
system. The mathematical model defining the input-output 
relationship of a physical system is termed the transfer 
function. 

By transforming the system excitation and system 
response to a complex domain, denoted by the subsidiary 
variable s, an exact expression for the transfer function 
is obtained. If the system has a single input and single 
output, the transfer function can be represented by means 
of a block diagram as shown in Figure 2. When expressed in 
equation form, 

X(s) - H(s)F (s) 2.12 

the transfer function can be defined as an operator. 
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x 


Transformed 

Input 

F(s) 


Transfer 

Function 


Hcs) 


Transformed 

Response 


X(S) 


Figure 2. Block Diagram for a Transfer Function 


In order to demonstrate the transfer function method 
for system characterization, consider the system defined in 
the previous section. If equation 2.3 is transformed to 
the Laplace domain, the following expression results: 

s 2 X(s) + 2Cw n aX(s) + w 2 X(s) - ^ F(a) 2.13 

The system transfer function is found by rearranging equa- 
tion 2.13 to the form 


H(a) m Mil 

H(a) F (s) 


1 

m 


2 2 
s + 2£w n s + w n 


2.14 


By substituting jw for s in equation 2.14 and rearranging, 
the following expression is obtained, 


H(w) - 


1/k 


2.15 


l-(o>/w n ) +j 2Cw/w n 
It is noteworthy that comparison of equations 2.15 and 
equation 2.11 prove them to be identical, the complex fre- 
quency response. 


It is evident that the transfer function method is 
another technique to evaluate the frequency response data 
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of a system. The transfer function of a complex system can 
be determined from experimental data. Thu®/ a system can 
be identified by data from a frequency response test, 

2.2 Modal. Analysi s 

An alternative method for characterizing a system's 
dynamic behavior is by defining each of the system's com- 
ponent modes of vibration. A discrete dynamic system has 
as many natural frequencies and modes of vibration as 
degrees of freedom. The general motion of the system can 
be described by the superposition of the modes of vibra- 
tion. 


Modal analysis is defined as the process of character- 
izing the dynamics of a system in terms of its modes of 
vibration. Each mode is described in terms of three modal 
parameters: an undamped natural frequency, a measure of 
energy dissipation, damping, and a characteristic deflec- 
tion shape, mode shape. A mode shape is a unique deforma- 
tion shape that the system would acquire if excited solely 
at the frequency associated with that mode. Determination 
of the three modal parameters for the modes of a system is 
the task of modal analysis. 

2.2.1 Analytical Modal Analysis 

In order to mathemtically interpret the modal parame- 
ters, a simple structural model is analyzed. Consider a 
linear, undamped system of n degrees of freedom represented 
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by the differential equations of motion in matrix form [4], 

[mni’(t)) + [k] (x(t) } - {0} 2.16 

If solutions to the equation are assumed to be of the form 

X| ( t) - u i e Xt 2.17 

and 

x t (t) - X 2 u i e Xt - X 2 X;L (t) 2.18 

then by substituting into the equations of motion and rear- 
ranging, the following equation results: 

[^[m] + W)(u) - {0) 2.19 

Equation 2.19 is a set of simultaneous algebraic equations 

2 

in u, the unknowns being the u's and X . For a nontrivial 

solution, {u} i* 0, the determinant of the coefficients must 

equal zero. The determinant of the coefficients is a poly- 
2 

nomial in X . The roots of the polynomial are called the 
eigenvalues,, If the eigenvalues are imaginary, then the 
modulus of the eigenvalues are equal to the natural fre- 
quencies of the system. 

The solution vector, {u}, corresponding to a particu- 
lar eigenvalue is termed an eigenvector. The eigenvectors 
are also referred to as modal vectors and represent physi- 
cally a deformation pattern of the structure for a natural 
frequency of vibration. Since equation 2.19 is homogene- 
ous, the solution for the eigenvectors is not unique. The 
vectors are unique only in the sense that the ratio between 
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any two elements in the vector is constant. Hence/ the 
shape of the natural modes is unique, but not the ampli- 
tude. The process of adjusting the elements of the natural 
mode to render the amplitudes unique is normalization, and 
the resulting vectors are referred to as normal modes. 

The natural modes satisfy a property called ortho- 
gonality. The statements of orthogonality with respect to 
the mass and stiffness matrices are as follows: 

{u>t[f0<u} 8 - (0} 2.20 

and 

{u)^tH](u} 8 - (0} 2.21 

where r and s specify eigenvectors of different eigenvalue 
solution. The proof of the orthogonality of the modal vec- 
to appears in the text and will not be demonstrated here 
[ 5 ]. 

The orthogonality property serves a useful purpose by 
allowing the equations of motion to be decoupled. For 

example , the modal matrix can be formed from the modal vec- 

tors in the form 


[u] - [{u} 1 {u} 2 ...{u} n ] 2.22 

The modal matrix functions to transform the mass and stiff- 
ness matrices to the generalized mass and stiffness 
matr ices , 

[u] t [m][u] - [M] 


2.23 
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and 

[u] t [k][u] - [K] 2.24 

The generalised mass and stiffness matrices are both diago- 
nal matrices due to the transformation. 

If a new coordinate system is defined by the transfor- 
mation 


{x} - [u] {q} 2.25 

then, substituting into equation 2.16, the equations of 
motion take the form 

Cm][u]{q‘} + C k] [u] { q } - {0} 2.26 

Pr e-mult iplying equation 2.26 by [u^ and substituting from 
equations 2.23 and 2.24 reduces the equations of motion to 

[M] {q " } + [K] { q } - {0} 2.27 

Since the general mass and stiffness matrices are diagonal, 
the coordinate transformation has completely decoupled the 
equations of motion. Each equation in equation 2.27 is 
equivalent to a single degree of freedom and is easily 
solved by the method of Section 2.1.1. 

The process of reducing a set of coupled system equa- 
tions to an uncoupled set by coordinate transformation and 
performing analysis in tne modal space from an nalytical 
sense i3 known as modal analysis. The same mathematical 
principles can be applied to systems with forcing functions 
or damping. When the damping is proportional to mass and 
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stiffness, the reaulting coordinate transformation reaulta 
in a aet of uncoupled force reaponae relationahipa exactly 
of the form of equation 2.3. 

2.2.2 Experimental Modal ^Analysis 

With the advent of the digital computer and tho intro- 
duction of an efficient Fourier integral algorithm in 1965, 
experimental modal analyaia haa evolved into a frequently 
uaed and important engineering analyaia tool. The theory 
governing atructural modal analyaia will be briefly 
reviewed aa background for application to acouatical modal 
analysis . 


Recalling the damped, 3 ingle degree of freedom ayatem 
from Section 2.1, it can be demonatrated that the transfer 
function, equation 2.14, can be rewritten to the form 

H(a) - — — 2,28 

(a-d) (s-d ) 

where d represents the pole of the transfer function and is 
given by 


d - <u> n ± I l-C 2 2.29 

and * is the complex conjugate operator [6]. Expanding 
equation 2.28 in partial fraction form results in 


H(s) - 


1/m 

(s-d) (s-d*) 


(S-d) 


* 

(s-d ) 


2.30 


The residues of the transfer function are defined by the 


constants and R 2 - 


The residues are directly related to 


13 


the amplitude of the response function. It is found by 
solving equation 2.31 that R 1 and are a complex conju- 
gate pair. Thus, equation 2.30 becomes 

H(s> ■ 75?37 + 2 ‘ 31 

The concepts of poles and residues associated with a 
transfer function are integral to experimental modal 
analysis and will be recalled later. 


An alternate method for determining the mode shapes 
from the eigenvalue problem of the previous section is by 
the use of an adjoint matrix [7]. If equation 2.19 can be 
reformulated to the expression given by 

[[k]' 1 [m]-X 2 [i]]{u) - {0} 2.32 

then an additional matrix i3 defined by 

[B(X)1 - [[k] -1 [m] - x 2 [lj] 2.33 

From the identities 

CB(X)][B(X)3" 1 - [I] 2.34 

and 


[B ( X ) ] 


2.35 


-i „ lean 

IB(X) j 

where [(C(X)J is the adjoint of matrix [(B(X)], the eigen- 


value problem for solution, r, takes the form 


If X 


2 

r 


CB(X r )][C(X r )] - |B(X t )|[I] 2.36 

is a root of the characteristic equation 2.32, then 
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| B ( X r ) |-0 2.37 

and equation 2.36 becomes 

[B(X r )J[C(X r )] - [0] 2,38 

Rewriting equation 2.38 using only the j-th column of 
[C(X )] yields the equation 

[B(X r )]{C(X r )} j - {0} 2.39 

Equation 2.39, like equation 2.32, represents a set of 
homogeneous equations in u which determine each value to 
within an arbitrary constant. Thus, vector {C(X r )}^ and 
{ u } are proportional and represent the mode shape associ- 
ated with the eigenvalue r determined from the characteris- 
tic equation. 

Similarly, consider the Laplace transform of equation 
2.16 with the addition of forcing functions 

[cm]8 2 + [3c] ] (X(8) ) - {F (3) } 2.40 

By substituting 

[B ( s ) ] = [[m]s 2 + [k]] 2.41 

into equation 2.40, the equations of motion become 

[B(s)]{X(a)} - { F ( s ) } 2.42 

where [B(s)] is referred to as the system matrix. Pre- 
multiplying equation 2.42 by [B(s)] 1 yields 

[B(s)] _1 {F(s)} - 'X(s)} 2.43 

By specifying the relationship between the transfer 
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function matrix and the system matrix, 

[H(8)] - [B(s)] _1 - 2.44 

the system responses are related to the system forcing 
functions through the transfer function matrix. Evaluating 
the adjoint matrix for the eigenvalues determined from the 
characteristic equation, it is apparent that the transfer 
function matrix contains both the eigenvalues and the 
eigenvectors to a proportionate constant. 

The method for uncoupling the equations of motion by 
the adjoint matrix when applied to systems with propor- 
tionate damping results in the same adjoint matrix as for 
the undamped system. The poles of the characteristic equa- 
tion are no longer the undamped natural frequencies, how- 
ever, and are given in the form of equation 2.29. 

The transfer function matrix is represented in the 
form [8] 


11<S) 

H 21 (8) 

H ml (3 > 

12< 6) 

« 22 (S) 

H n,2 (9 > 

l m (s) 

H 2m < 8) 

%, (s) 


where m is the number of degrees-of-f reedom. Substituting 
equation 2.45 into equation 2.43 and multiplying through 
yields 


H 11 ( 8 )F 1 (8) + --' + H lm (3)F m (s) * X l (s) 


2.46 
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H ,(s)F,(s) + ...+ H _ ( a ) F ( 8 ) * X m (s) 
ml ''l' 7 mm' ' m ' m' ' 

If in equations 2.46 all forcing functions are set to zero 
except F 1 (s), then the individual transfer functions become 


x i(a) 

H ll< 9 > " iT7i7 


.47 


H 21 (s) - 


X 2 (s) 

F^TiT 


X ( s ) 

H ml (s ^ “ F^TsT 

Thus, the first column of the transfer function matrix 
relates the response at all m degrees-of -freedom to an 
excitation at degree-of-f reedom 1. It is noteworthy that 
all the individual transfer functions have the same chara- 
teristic equation. 


Each of the transfer functions in equation 2.45 can be 
represented as the sum of the transfer functions for single 
degree-of-f reedom systems in the form of equation 2.31, 


H ab (3) - 

r-1 


R r (a,b) R r (a,b) 
+ i 


(s-d r ) 


(s-d ) 
r 


2.48 


where H &b (s) relates the response at degree-of-f reedom a to 
an excitation at degree-of-f reedom b, and n is the number 
of roots to the characteristic equation. Thus, any of the 
transfer functions in equation 2.47 can be used to find the 
complex poles of equation 2.48. The transfer function 
matrix is now expressed in the form 
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C H ( s ) ] 


n 

E 

r-1 


CV 


[R*] 


(3-d ) 
r 


2.49 


(s-d r ) 

Comparing equation 2.49 and equation 2.44, it is apparent 
that the residue matrix is equivalent to the adjoint 
matrix. Prom the definition of the system adjoint matrix, 
the residue matrix, [R r ], evaluated at the natural frequen- 
cies is proportional to the mode shapes by the relationship 


[R r ] - « r {u r }(u t } t 2.50 

Thus, each column of the residue matrix evaluated at the 
natural frequencies is proportional to the eigenvectors. 

The three modal parameters , <*> r , C r , and u r , are found from 
a single row or column of the transfer function matrix. 


2.3 Acoustical Modal Analysis 


2.3.1 Literature Survey 

The analytical and experimental study of reactive 
acoustic fields and resonant modes has progressed, and 
record exists in the literature. A survey of recent 
developments will contribute to defining the direction of 
this research. 


In some early work, Gladwell applied variational prin- 
ciples to formulate certain acoustical and structural con- 
figurations [9]. Later, Craggs considered variational 
statements in the context of finite element representations 
for the coupled structural-acoustic problem [10]. A 
comprehensive structural-acoustic theoretical model was 
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developed for interior sound fields using modal procedures 
by Dowell, Gorman, and Smith [11]. The sound fields which 
are created by arbitrary wall motions were modelled using 
Green's Theorem. Petyt, Lea, and Koopmann developed a fin- 
ite element method for determining the acoustic modes of 
irregular cavities [12]. In the same time period, Wolf, 
Nefske, and Howell combined Dowell's modal procedures with 
both structural and acoustic finite element analysis to 
study automobile passenger compartments [13j. 

The finite element method lends itself readily to the 
acoustic problem, and many engineering problems are being 
analyzed by this method. Finite element modeling does, 
however, suffer two major disadvantages. Most acoustic 
finite element computer programs are very large in size, 
and thus require large computers with ample memory in which 
to operate. In order to obtain, the required accuracy, 
models containing several thousand degrees-of-f reedom are 
not uncommon. A second disadvantage of finite element 
modeling i3 the error due to the use of an inadequate 
number of elements or unrealistic boundary conditions. 

It is particularly true in structural analysis that 
frequency response can be closely approximated using a 
truncated set of the total set of modes (usually the lower 
frequency modes) . By using the procedures of equations 
2.24 through 2.28, a significantly smaller approximate 
model can be developed. Thus, modal analysis where 
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acoustic and/or structural modes are found either experi- 
mentally or analytically and coupled to build an overall 
structural-acoustic modal model are very popular. 

The technology for experimental modal analysis of 
acoustic fields is lagging both finite element modeling and 
experimental structural modal analysis. The majority of 
the research performed has been for the one-dimensional 
problem. Ibrahim and Mikulcik devised an experimental 
technique by which standing wave parameters in gas piping 
systems could be identified in the time domain [14]. Two 
pressure taps were made in an impedance tube. The tube was 
closed at one end and excited by a sweeping sine wave at 
the opposite end. The modal parameters were computed from 
the decay rates and amplitude ratios between the micro- 
phones. Experimental results agreed well with the theoret- 
ical results obtained from the undamped, plane wave equa- 
tion. 


A technique for experimental modal analysis of cavi- 
ties involving Fourier analysis was proposed by Smith [15]. 
Evaluation of a one-dimensional tube with a small volume 
velocity source as the input showed that the response sound 
pressure is in-phase with the source volume velocity only 
at a cavity resonant frequency. Smith determined the 
acoustic pressure mode shapes by measuring transfer func- 
tions for several microphone locations and plotting the 
magnitude and phase at each resonant frequency. 
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In the results of a recent study, Mehl contended that 
background noise in pressure mode shape measurements was 
significant [16]. Mehl described numerical techniques for 
eliminating the effects of background noise in precise 
measurements of normal mode frequencies and damping. 
Resonant standing wave measurements were fitted to theoret- 
ical resonance formulas. 

Neiter and Singh describe an acoustic modal analysis 
experiment employing existing software for structural modal 
analysis [17]. They describe the inherent characteristics 
of the acoustic system in terms of an acoustic impedance 
matrix, Z, by the relationship, 

(P(s) } - [Z(s)]{Q(s)} 2.51 

where Q is the complex volume velocity source and p is the 
complex pressure response.. By curve-fitting the measured 
impedance matrix to a frequency response function similar 
to equation 2.48, the modal parameters were extracted. The 
experimental results agreed well with a theoretical model 
possessing no damping. 

In a report for research performed at Ray W. Herrick 
Laboratories, Byrne describes a procedure for experimen- 
tally determining, in terms of particle motions, the shapes 
of the low-order acoustic modes in three-dimensional enclo- 
sures [18]. Transfer functions were measured relating 
acoustic pressure response at an array of points in the 


enclosure to a volume velocity source. The pressure mode 
shapes were described discretely by the peak amplitudes 
extracted from 1 the transfer functions. A differentiable 
function was curve-fit to the discrete pressure amplitudes. 
The gradients of each approximate pressure mode function 
are determined and related to particle motion by the invis- 
cid force equation, 

P 0 |a - -Vp 2.52 

Thus, Byrne was able to describe the acoustic mode shapes 
in terms of particle motion. 

2.3.2 Rg.B ^ .ai.g±L-gJjLg-<? t j Qh 

The primary objective for the research presented here 
is the development of a procedure by which the three- 
dimensional modal response of an acoustic cavity is deter- 
mined experimentally. It is advantageous to describe the 
acoustic mode shapes in terms of particle motion for exci- 
tation source identification. In light of the research 
survey presented in the previous section, the existing 
technology founded on the theory of modal analysis will be 
applied to the objectives. In this section, the research 
direction for application to the objectives is outlined. 

A method for extracting the modal parameters of a 
mechanical system from the complex frequency response 
matrix is outlined in the first two sections. The system 
displacement response is related to the forcing functions 
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by the system matrix as given in equation 2.42. The analo- 
gous acoustical relationship for determining the complex 
pressure response due to a volume velocity source is 
expressed by equation 2.51. By assuming a lightly damped 
system, decoupled modes, and no background noise, the modal 
vectors are obtained directly from the residues of the 
measured frequency response functions. 

The frequency response method assumes a lumped parame- 
ter system with a modal vector element for each degree of 
freedom. The mode shapes associated with the natural fre- 
quencies of a continuous system, however, are waveforms. 

For an acoustic system the waveforms correspond to standing 
waves. When the frequency response method is applied to 
continuous systems, the modal vectors become a discrete 
representation of the waveform deflection shapes. A func- 
tion which approximates the mode shapes is obtained by per- 
forming a curve-fit procedure to the modal vectors. In 
acoustics, if the approximate pressure mode shape function 
is differentiable in spa/::e, then the mode shapes can be 
described in terms of particle motion by the inviscid force 
equation, equation 2.52. 

In order to satisfy the objective for the development 
of an acoustic modal analysis procedure, the directions for 
research are outlined. The structural modal analysis sys- 
tems commonly in use are all-inclusive mini-computer based 
systems. An additional research objective to be considered 
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is to tailor the acoustic modal analysis procedure includ- 
ing the curve-fit technique and graphic display to a mini- 
computer based spectral analysis system. The objective is 
accomplished by minimizing computation and storage space. 
Thus, the development of the acoustic modal analysis pro- 
cedure will be performed on an existing spectral analysis 
system. 


2.3.3 Theoretical Model 

A theoretical model provides a reference for develop- 
ing and verifying the acoustic modal analysis procedure. 

The model is obtained from solutions to the well-known 
acoustic wave equation, expressed in vector notation in 
terms of acoustic pressure as follows [19]: 

2 1 8^p 

V p - — t — £ 2 . 53 

c at 

Equation 2.53 is the linearized, lossless wave equation for 
the propagation of sound in fluids. The general harmonic 
solution to the wave equation is of particular interest. 

If a solution in three-dimensional Cartesian coordinates of 
the form, 

p(x,y, z , t) - F(x,y,z)e Xt 2.54 

is substituted into equation 2.53, then by a separation of 
variables procedure the solution is found to be [20] 

p(x,y,z,t) - X(x)Y(y)Z(z)e Xt 2.55 


where 


- - - u 


X(x) - A cos k x + B sin k x 
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2.56 * 


X(y) - c cos k y y + D sin k^y 


Z(z) - E cos k z + F sin k z 
2 2 

In order to model, an acoustic cavity, it is necessary 
to specify the domain boundary conditions. Acoustic cavity 
t e#onances are indicative of boundaries of relatively high 
acoustic impendance. For the free response problem, the 
acoustic cavity is assumed to be of dimensions 
L x by Ly by L z with orthogonal boundaries of infinite 
acoustic impendance located on the major axes and at 
L x' L y' and L z' fundamental acoustic principles at the 
boundaries , 


n*u - 0 
n'vp - o 

Thus, the boundary conditions bee me 


(Ifi) - 0 

v dx'x«0,L x 

(•§£) . o 

WV-OrLy 


2.57 

2.58 


2.59 


(4») - o 

k 3z'z=0,L z 

When the boundary conditions are applied to the general 
harmonic solution, the following expession results 


± 2 * cos MX co8 ffi e ^ lmnt 
x y z 


’ p lmn coa 


z . so 
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The symbol X represents the eigenvalues, and the coeffi- 
cient ' p lmn' 18 com P lex - The abeolute value of imaginary 
eigenvalues are equal to the natural frequencies and are 
given by the relationship 


,, wr/ 1.2 , . m. 2 , . n. 2,1/2 

l X l - u lmn - c7r t (l - ) + <iT> + <L~> 3 

x y z 

where 1, m, and n are the respective mode numbers. 
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The solutions to the theoretical model, given by equa- 
tion 2.60, represent three-dimensional standing waves in a 
rectangular cavity. In order to characterize the system by 
modal analysis techniques, the response is expressed as the 
sum of the individual .pressure mode shapes in the form 

00 00 00 

p(x,y,z,t) - Z Z Z p lmn* cos iT* 

1-0 m-0 n-0 lmn x 


cos 


mfry 

L y 


cos 



^lmn^ 


2.62 


By determining the gradient of the individual pressure 
mode functions with respect to the three coordinate direc- 
tions, the particle acceleration associated with the 
respective pressure mode is found via the inviscid force 
equation, equation 2.52. The gradient operator in Carte- 
sian coordinates is given by 


v “ a x + 3y + 3z 2,63 

The partial derivative of the pressure mode function, equa- 


tion 2.60, with respect to the x coordinate is found to be 
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6x 



lmn 


cos 


mix 

L y 


COS 


n TT'Z 
L_ " 


sin 


It rx 


lmn 


2.64 


The partial derivatives of the pressure mode function with 
respect to the y and z coordinates are found similarly. 
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3. CURVE-PIT TO THREE-DIMENS lONAL DISCRETE DATA 

From the frequency response method presented in the previ- 
ous section, the modal vectors are obtained directly from 
the residues of the frequency response functions. The 
modal vectors resulting from the acoustic modal analysis 
procedure established in Section 2.3.2 discretely represent 
three-dimensional pressure standing waves. Functions which 
approximate the acoustic mode shapes are obtained by per- 
forming a curve-fit in space to the modal vectors. In this 
section a procedure for fitting a function to three- 
dimensional discrete data is developed. 

An objective for the acoustic modal analysis procedure 
is to describe the acoustic mode shapes in terms of parti- 
cle motion. From the previously cited inviscid force equa- 
tion, it was established that the acceleration of a gas 
particle is given by the negative of the gradient of the 
pressure field at the location of the particle. Therefore, 
the particle acceleration mode shapes are obtained if the 
acoustic pressure mode shapes are described by a differen- 
tiable function. An objective for the curve-fit procedure 
is to derive an approximate pressure mode function which is 
differentiable with respect to the spatial coordinates. 

The character of the acoustic mode shapes in an enclo- 
sure is directly dependent on the dimensional shape of the 
cavity and the surface boundary conditions. Thus, the most 
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that can be assumed concerning the pressure function 
describing an acoustic mode in an enclosure is that the 
shape is sinusoidal in space, due to the physical nature of 
sound. In addition, each data point in the modal vector is 
subject to experimental error which, for acoustic measure- 
ment, may be relatively large in magnitude. It is inap- 
propriate to match an interpolating function exactly at the 
data points since no implicit assumption can be made con- 
cerning the accuracy of the data. The pressure mode func- 
tion, therefore, must be assumed arbitrary. Thus, the 
objective is to find the 'best' curve which represents data 
that are subject to error. The criterion for goodness of 
fit is to some degree arbitrary, although the least squares 
criterion is most common and will be applied to the acous- 
tic problem. The least squares procedure offers the bene- 
fit of a polynomial interpolating function which is numeri- 
cally differentiable. 

3 . 1 Least Squares Procedure 

In order to define a three-dimensional acoustic 
domain, the Cartesian axes system is specified for this 
study. Thus, the function describing the acoustic pressure 

1. U 

field associated with the r acoustic mode in an enclosure 
is p (x,y,z). The procedure objective is to derive an 
approximate function to P r (x,y,z) which is differentiable 
with respect to the spatial coordinates. If the approxi- 
mating function is specified to be h (x,y,z), then thf* 
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following relationship holds 

P r (x,y,z) - h E (x f y,z) + e(x,y,z) 3.1 

The term e(x,y,z) represents error associated with the fit. 


The least squares procedure developed by Gauss can be 
stated verbally as, 'the sum of the squares of the residu- 
als of the differential equations should be minimum at the 
correct solution' [21]. The least squares error norm 
expressed mathematically for a function in three dimensions 
is as follows [22] 



where from equation 3.1 


J//e(x r y r z) 


2 


dxdydz 


3.2 


e(x,y,z) - p r (x,y,z) - h r (x,y,z) 3.3 

It is now necessary to expand the approximating function 
into a series of interpolating functions and constants 

h r (x,y,z) - c 1 <P 1 (x,y,z) +...+ c i <P i (x r y,z) 3.4 

where i represents the order of the function. The error 
term is now given by 


e(x r y r z) = p r (x,y,z) - c 1 «t> 1 (x,y,z) -... 

- c i 0 i (x,y r z) 3.5 

and the least squares error norm becomes 

2 

Mb " •f-^Pr “ c l*l c i <p i^ dxdydz 3.6 

In order to complete the least squares procedure and 


minimize the error associated with the curve fit, the 


partial ‘derivative of equation 3.6 with respect to each 
constant is determined and set to zero 
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0 “ SSS 2 [p r - c^ 


c^] (-^dxdydz 3.7 


3 e 


|/3 Ci - 0 - SIS 2 [p r - c-^ -...- c i <P i ] (-1» i )dxdydz 
Equation 3.7 can be reduced to the form 

2 

SSSp t ^dxdydz * c^/J*/<b^dxdydx + ,..+ c^/J* /(b-^^dxdydz 3.8 


//Jp r ^idxdydz = c^/JJV^dxdydz +...+ c,. J//<l> i dxdydz 
Equations 3.8 are rearranged to matrix form as 


;;;*2 axdydz 

///(b^^^dxdydz 


C 1 


///Pr^ldxdydz 

X 

X 

< 

X 

) •* < 

x ) 

J*J*J*<b ^ ^ d xdy d z 

_ 

J*J* j*<i>^dxdydz 


c , 
1 


SSSp Y ^.dxdydz 

l J 


>3.9 


By solving equations 3.9 for the constants and substituting 
those valves into equation 3.4, the function, h r (x,y,z) , 
has been determined which 'best' approximates the curve of 
the discrete data with respect to the least squares error 
norm. 


3 . 2 Legendre Polynomials 

A critical element in the least squares procedure is 
specifying the interpolating function to approximate the 
acoustic pressure mode function. In order to select the 
type of function most appropriate to the application, it is 
necessary to consider some important criteria: 


H 


( 
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a. The accuracy of the approximation is directly 
dependent on how well the polynomial function can 
approximate acoustic data as evaluated by some 
error norm. It is noteworthy that acoustic 
phenomena are sinusoidal in nature. 

b. In order to solve the matrix equations for the 
constants of the polynomials, the polynomial 
function must be integrated with respect to the 
spatial coordinates. The ease at which the 
integration is performed becomes a key element of 
the procedure. 

c. Recalling the objectives from the introduction, 
the acoustic modal analysis procedure should be 
tailored to mini-computers. Therefore, computa- 
tion efficiency is greatly desired. 

After a review of the available polynomial shape func- 
tions conducive to the least squares procedure, Legendre 
polynomials are selected for approximating the acoustic 
data. With reference to the aforementioned criteria, the 
selection was made on the following premise: 

a. All complete interpolating polynomial families of 

fch 

sequential n order will approximate a 
sinusoidal function with equivalent accuracy. 
Thus, no accuracy is sacrificed by the use of 
Legendre polynomials over another interpolating 
polynomial . 


32 


b. Legendre polynomials are separable with respect 
to the three spatial coordinates. Therefore, the 
polynomial elements are integrated with respect 
to the associated spatial coordinates rather than 
being integrated over a unit volume element. 


c. Legendre polynomials are orthonormal in the range 
from -1 to 1. The significance of this charac- 
teristic is that if all functions of interest to 
the least squares procedure are defined on the 
region -1 to 1, then all off-diagonal terms of 
the left-side matrix become zero. Thus, a con- 
siderable reduction in computation time has been 
attained. 


The Legendre polynomials are 


S Q (u) 

= 1 

S 3 (u) 

V u) 

= u 

S 4 (u) 

s 2 (u) 

- 1/2 ( 3u 2 -1) 

V u) 


„ , v 2n-l 



S n u ■ n 

uS n-l ( 


listed as follows [23].* 
1/2 ( 5 u 3 -3u) 

1/8(35u 4 ~30u 2 +3) 3.10 

1/8(63u 5 -70u 3 +15u) 


In order to approximate a function of three- 
dimensions, it is necessary to create shape functions which 
are a product of three Legendre polynomials in the form 


~ s i x < x > s my (y) s nz (z > 


3.11 
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where l,m,n represent the polynomial order in the respec- 
tive directions. The maximum order required for the Legen- 
dre polynomials in the three directions is an important 
consideration and will he addressed in a subsequent sec- 
tion. The order of the shape function, i, is dependent on 
the order of the three Legendre polynomials and becomes a 
matter of bookkeeping when programmed in FORTRAN. 


Recalling that Legendre polymonials are orthonormal in 
the region -1 to 1, the three dimensional acoustic domain 
must be mapped to a 2 x 2 x 2 cubic element in the C, 7?, £ 
domain. As a result of the Legendre polynomials being 
separable with respect to the three coordinate directions, 
and since the acoustic domain is defined on an orthogonal 
axes system, the map reduces to merely a scaling in each of 
the three directions. The scale factors only become impor- 
tant when determining the gradient of the approximation 
function. The matrix equation, equation 3.9, is now 
transformed to 


/ 




< 


1 1 1 2 

///<*>, (C ,7?,OdCd7?d€ 
-1-1-1 


0 


0 


1 1 1 , 

SIS 4>f(C ,7?,£)dCd7?d£ 
- 1 - 1-1 1 


111 

S S S p r 4> (C,7?,£)d<d7?d* 
- 1 - 1-1 r 1 

< x > 

111 

SIS p (C,7?,«)dCd7?dS 
- 1 - 1-1 r 1 


c . 

1 


3.12 
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Computation can be further reduced by modifying the 
left aide matrix in order to create the identity matrix. 

By obtaining the identity matrix from the left side of the 
equation, the constants, c i , are solved for directly. An 
analysis of Legendre polynomials reveals that if the square 
of the polynomial is integrated exactly over the limits -1 
to 1, a number results with a value directly associated 
with the order of the polynomial 

1 2 

J S n (u) du - 2/2n+l 3.13 

-1 n 

Thus, a series of constants can be defined to compliment 
the respective shape functions 

1 1 1 2 111 2 

= -i^i tS a« (<)S M (,,,S c « (i)1 ‘ Kd,W * 

- (2/2a+l) (2/2b+l) (2/2c+l) - P t 3.14 

When both sides of the matrix equation 3.12 are divided by 
the constants, p^, the left side matrix becomes the iden- 
tity matrix, and the right side vector acquires the 
reciprocal of the normalizing constant 


1 0 

' C 1 


/111 \ 
1/Pi/ S I p r <M£,7?,*)d<:d7?d£ 

x xL 

X 

> * ( 

-1-1-1 \ 
1 1 1 1 ! 

1 ° *1 

c i 


\i/pJ S S p*At,'nfZ)acavdz 
\ 1 ~'i—i—i r 1 J 


3 . 3 Numerical Integration 


The matrix equation derived by the least squares pro- 
cedure has now reduced to a series of integral equations: 
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r ; 


111 

c, - l/P-J / J* p r <P 1 (C,7?,€)dCd7?di 3.16 

- 1 - 1-1 

111 

c t - l/PJ I I P-Pj ( C ,v , Oc3Cd7?di 

The constants that yield the best curve fit with respect to 
the least squares norm can be solved for directly by per- 
forming the required integration. In order to adapt the 
curve fitting procedure to a computer, the computing struc- 
ture requires that the integration be performed numeri- 
cally. Several numerical integration methods are in 
existence. The criterion for the method of integration is 
again accuracy and efficiency of computation. 

Consideration of the existing numerical methods of 
integration results in the selection of the Newton-Coates 
numerical integration formulas. The Newton-Coates integra- 
tion equations have the general form [24] 

x n 

I f(x)dx - K n h(w 0n f 0 + vt ln f 1 +...+ w nn f n ) 3.17 

x 0 

where 

K n - the n-point scheme constant 
h - the distance between points 
f = the value of the function at the points 
w * the weighting constants of the function valves 
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The accuracy of the integration can be maximized by apply- 
ing the point scheme equation that coincides exactly with 
the number of recorded values of the function being 
integrated. Alternatively, a combination of Simpson's 1/3 
rule, 3/8 rule, and the trapezoidal rule can be applied for 
integration over any number of points. 

By applying the Newton-Coates formula, equations 3.16 
become 

nC M) n£ 

C 1 - l/p x t t l P rabc P 1 (< a -V« c )w a w b w e 3.18 
a-i b-1 C“1 

nC nv n£ 

c i ‘ VPi a fl b -l c f/^abcV i a'V { c> w a w b w c 
The numerical integration can now be performed, and the 
constants of the approximation polynomial function are 
solved for directly. The structure of the mathematical 
procedure lends itself readily to FORTRAN or BASIC program- 
ming on a mini-computer. 

3 . 4 Error Analysis 

An important consideration for the development and 
evaluation of the modal analysis procedure is the accuracy 
with which the computed pressure mode function approximates 
the discretely measured pressure mode shapes. Recalling 
that the pressure mode shapes are sinusoidal in space, 
three questions concerning accuracy of fit should be 


addressed : 
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a. What is the minimum polynomial order that suffi- 
ciently approximates a sine function of a given 
number of wavelengths? 

b. How many data points per wavelength are required 
to adequately describe a sine function? 

c. What is the error associated with a numerical 
integration? 

An analysis of the curve-fitting procedure with respect to 
the three points of question will allow for the design of 
an optimal configuration. 


In order to evaluate the accuracy of the numerical 
procedure, a known function is approximated and an error 
criterion is established. After consideration of the vari- 
ous error norms, a standard deviation error criterion was 
chosen for judging accuracy. This error norm was chosen due 
to the discrete nature of the data and ease of evaluation. 
The error norm expressed in equation form is [25] 


e 


n 

£ 

i=l 


(P r (i)~h r (i )) 2 


1/2 


3 * 19 


3.4.1 Polynomial Order 

In order to evaluate error as a function of polynomial 
order, the error due to a polynomial approximation must be 
segregated from error due to other variables. By assuming 
a function, p r , and performing an exact integration of 
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equation 3,16, the approximate function, h r , is computed 
independent of numerical integration err m or number of 
data points. 

Recalling from the theoretical model, the pressure 
mode shape for the ideal case in one dimension will be of 
the form 

p r (x) - cos 1 ttx/L x 3,20 

Equation 3.20 represents the one-dimensional pressure mode 
shapes defined for the domain 0 < x < L v . In order to com- 
pute the coefficients of the approximation function, L is 

X 

set to 2 and Equation 3.16 is modified such that the limits 
of integration are 0 and 2, By substituting (u-1) for u in 
the Legendre polynomials, equation 3.10, the limits of 
integration can be adjusted appropriately. With substitu- 
tion of equation 3.20 and the modified Legendre polynomials 
into equation 3.16 and exact integration, the approximation 
function can be found. 

The approximation function was computed for the first 
four modes, 1 * 1,2, 3, 4, with approximation functions of up 
to and including eighth-order polynomials. The error asso- 
ciated with each case is evaluated for eleven sample points 
via equation 3.19, and the results are presented in Table 
1. The results show that in the absence of integration 
error, the approximation improves with increasing polyno- 
mial order. 
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3.4.2 Numerical Integration 

The coefficients for the approximation function result 
from the numerical integration in equation 3.18. In order 
to perform the computation, the pressure mode function for 
one-dimension, equation 3.20, is described in terms of 
discrete values distributed evenly along the x-axis. It is 
intuitively obvious that the approximation error will 
approach zero as the increment between points decreases and 
the integration approaches an exact integration as dis- 
cussed in the last section. A study of the error as a 
function of the number of integration points for an ideal- 
ized case, however, illustrates the sensitivity of the 
approximation technique to integration error. Such error 
may lend insight into the design of later experiments. 
Therefore, rather than assume to minimize increment between 
points which increases data collection time, computation, 
and data storage space, the number of data points is a 
variable in the analysis. 

An inverse relationship between error and polynomial 
order was proven for an exact integration approximation. 

The polynomial order for the curve-fitting procedure must 
be limited, however, to the benefit of computation time. 

In addition, a polynomial approximation of a sine wave may 
behave peculiarly for low order polynomials due to numeri- 
cal integration error. The polynomial order, therefore, is 
an important variable in the error analysis. 


i* ' 
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The coefficients of the approximation functions 
corresponding to each combination of data points, polyno- 
mial order, and wave number is computed via equation 3.18. 
The error associated with each approximation is determined 
for eleven points in equation 3.19. The results are 
recorded and tabulated in Tables 2, 3, 4, and 5. 

The accuracy of the approximations are, in general, 
satisfactory. Certain trends in the error tables, however, 
are worth noting. The error is unacceptable for the cases 
of the polynomial order exceeding the number of integration 
points. A minimum of eleven integration points is required 
to adequately approximate a cos 1 1/2 wave (3rd mode). In 

4* u. 

addition, in order to approximate a cos 2 wave (4 mode) , 
greater than eleven integration points are required. 

Further evaluation of the error tables is performod for the 
design of the experiments. 

3 . 5 Approximation Function Gradient 

Resulting from the curve-fit procedure are polynomial 
functions which approximate the pressure fields associated 
with the acoustic modes in enclosures. From the inviscid 
force equation, equation 2.52, the acceleration of a gas 
particle is given by the negative of the gradient of the 
pressure field at the location of the particle. Thus, in 
order to describe the acoustic mode shapes in terms of par- 
ticle motion, the gradient of the approximate pressure 
function must be determined. 



Polynomial Order 


4 th 

5 th 

B 

7 th 


1.312 

1.312 

2.857 

2.857 

5.083 

.460 

.460 

2.428 

2.428 

3.629 

.122 

. 122 

1.078 

1.078 

2.866 

.068 

.068 

.535 

.535 


.004 

.004 
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.0038 
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or for a Full Cosine Vtove 


Polynomial Order 


4th 

5th 

6th 

7 th 

8th 

1.157 

1.425 

1.425 

4.687 

4.687 

.403 

1.511 

1.511 

2.480 

2.480 

.276 

.704 

.704 

1.8 

1.8 

.255 

.260 

.260 

1.35 

1.35 

.295 

.336 

.336 

1.047 

1.047 

.290 

.201 

.201 

.472 

.472 

.282 

.032 

.032 

.510 

.510 

.282 

.025 

.025 

.304 

.304 

.282 

.035 

.035 

.076 
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Irror for 1 1/2 Cosine Wave 


Polynomial Order 


4 th 

5th 

6 th 

7th 

8 th 

1.081 

1.081 

2.874 

2.874 

5.078 

1.433 

1.433 

1.433 

1.433 

3.998 

.765 

.765 

1.33 

1.33 

1.58 

.480 

.480 

.870 

.870 

1.913 

.436 

.436 

.757 

.757 

1.330 

.395 

.395 

.261 

.261 

1.116 

.515 

.515 

.780 

.780 

1.502 

.489 

.489 

.466 

. 466 

.604 

.447 

.447 

.145 

.145 
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Number of 


Polynomial Order 


Data Points 

19 

2nd 

3rd 

4th 

5 th 

6 th 

7 th 

8 th 

3 

1.225 

1.225 

1.225 

1.225 

2.646 

2.646 

3.338 

3.338 

4 

.801 

.801 

1.007 

1.007 

1.421 

1.421 

2.527 

2.527 

5 

.938 

.938 

.848 

.848 

1.604 

1.604 

1.042 

1.042 

6 

.820 

.820 

.755 

.755 

1.088 

1.088 

1.401 

1.401 

7 

.770 

.770 

.913 

.913 

1.002 

1.002 

.974 

.974 

8 

.769 

.769 

.778 

.778 

.617 

.617 

.731 

.731 

. 

9 

.778 

.778 

.705 

.705 

i 

.733 

.733 

1.252 

' 

1.252 

10 

.777 

.777 

.702 

.702 

.529 

.529 

.437 

.437 

11 

.774 

.774 

.717 

.717 

.851 

, 

.851 

. 

1.763 

1.763 
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The approximate pressure function results from substi- 
tuting the constant values obtained' in equation 3.18 into 
the following equation: 


h r (C,7?,*) - c 1 <P 1 ( C , i? f O + ...+ c.<P.(C,7?,0 3.21 

It is noteworthy that the approximate pressure function is 
defined in the C , V , i domain. In order to obtain the par- 
ticle acceleration for the original acoustic domain, 
defined in x, y, and z, the gradient of equation 3.21 with 
respect to x, y, and z must be computed. The partial 
derivative of equation 3.21 with respect to the x coordi- 
nate is found by the application of the chain rule as fol- 
lows 


ah r (4,7/,Q 

3x 


a h a h . an 

r £jC , r dji r 3J_ 

dC 3x 37? 3x 3£ 8x 


3.22 


Since both the x, y, z domain axes and the C , V , i domain 
axes are orthogonal, x is not a function of 7? or £ . Thus, 


3 * 


an 


and ^ are zero. Furthermore, the term is merely the 
3x 3x 3x 


constant ratio of the domain length along the C axes, 2, to 
the length of the acoustic domain in the x direction, L^. 
The same properties hold true for the partial derivatives 
of equation 3.21 with respect to the y and z coordinates. 


The variable terms resulting from the application of 
the chain rule to equation 3.21 represent the gradient of 
the approximate pressure function with respect to the C ; t? , 
and i coordinates. In order to compute the gradient, the 


derivatives of equations 3.10, the Legendre polynomials, 
with respect to u are determined. By substituting equation 
3.11 and equation 3.21 into equation 3.22, an expression 
for the partial derivative of the approximate pressure 
function results 


ah r <c,i>,€) 

3x 


[ 1 a s 1£ (o 

J 


' 2/ V 


3.23 


where 1, m, and n represent the Legendre polynomial orders 
associated with the shape function order, j. Expressions 
for the partial derivatives of equation 3.21 with respect 
to the y and z coordintes are found similarly to be 


and 


3h r (C,7?,€) 

dy 


3h r (C 

3z 


i a s _(t?) 

“ c j 8 K (0- #— 

3=1 J 


[2/L y ] 


3.24 


Lj f 1 C j S l^ <C)S mr 7 (7?) 


as n ,(0 


as 


[2/L ] 3.25 

Z 


Substitution of equations 3.23, 3.24, and 3.25 into the 
inviscid force equation yields an expression for the acous- 
tic mode shapes in terms of particle acceleration. The 
numerical operations in the three equations are conducive 
to FORTRAN or BASIC programming. 
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4. CONCLUSIONS 

The primary objective for the research presented here, as 
stated in the introduction, is the development of a pro- 
cedure for experimentally determining the modal parameters 
associated with resonance conditions in the cavities of 
transportation vechicles. It is desirable for noise con- 
trol diagnostics to describe the mode shapes in terms of 
particle motion. The acoustic modal analysis procedure 
should be adapted to function on existing computer based 
spectral analysis systems. 

In general, the objectives for the research are satis- 
fied. A procedure is developed whereby the acoustic mode 
shapes associated with the cavity resonances in an enclo- 
sure are determined experimentally and described in terms 
of particle motion. While not currently implemented on a 
mini-computer system, the modal analysis procedure includ- 
ing the curve-fitting program is constructed for use on 
small computer systems. 

The curve-fitting program performed satisfactorily for 
the modal analysis experiments by approximating the pres- 
sure mode functions discretely described by the modal vec- 
tors. Several improvements to the curve-fitting procedure 
are suggested, however. The procedure is limited due to 
the error resulting from fitting to a relatively few 
points. The program should be expanded to accomodate a 


larger number of data points . Furthermore, excessive error 
can result from extrapolating a polynomial function approx- 
imation of a sinusoid beyond the measurement domain. 

Either the acoustic pressure must be measured to the 
enclosing surfaces, or a shape function which approximates 
the harmonic acoustic pressure function beyond the measure- 
ment field should be introduced in the curve-fitting pro- 
gram. Finally, the requirement for equally spaced, orthog- 
onal input data to the curve-fitting program hampers meas- 
urements in irregular shaped domains. 

The task of data acquisition for the necessary number 
of points to adequately describe the acoustic modes in a 
typical transportation vehicle is long and tedious. 

Improved data collection techniques are a necessity. Two 
recommendations are for a microphone array attached to a 
multiple input spectrum analyzer or a stepper motor con- 
trolled point collection device. 
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